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We examine the effect of perturbative string loops on the cosmological pre-big-bang 
evolution. We study loop corrections derived from heterotic string theory compacti- 
^| fled on a Z N orbifold and we consider the effect of the all- order loop corrections to 
the Kahler potential and of the corrections to gravitational couplings, including both 
threshold corrections and corrections due to the mixed Kahler-gravitational anomaly. 
We find that string loops can drive the evolution into the region of the parameter space 
where a graceful exit is in principle possible, and we find solutions that, in the string 
frame, connect smoothly the superinflationary pre-big-bang evolution to a phase where 
the curvature and the derivative of the dilaton are decreasing. We also find that at a 
critical coupling the loop corrections to the Kahler potential induce a ghost-like insta- 
bility, i.e. the kinetic term of the dilaton vanishes. This is similar to what happens in 
Seiberg-Witten theory and signals the transition to a new regime where the light modes 
in the effective action are different and are related to the original ones by S-duality. In 
a string context, this means that we enter a D-brane dominated phase. 
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1 Introduction 



Pre-big-bang cosmology has been initially developed using the lowest-order effective 
action of the bosonic string This allowed to understand the basic features of 
this cosmological model: the Universe starts at weak coupling and low curvature, 
follows a superinflationary evolution, and enters a large curvature phase. As long as 
we are in the perturbative regime, a description with the lowest order effective action 
is adequate, and can be used to discuss the cosmological evolution [|IJ, the generality 
of the initial conditions @, and even to develop some interesting phenomenological 
consequences concerning the generation of primordial gravitational, axionic, dilatonic 
and electromagnetic backgrounds ||. In particular, the low frequency part of these 
spectra is only sensitive to the low curvature part of the evolution. 

At lowest order, the cosmological evolution unavoidably reaches large curvatures and 
strong coupling and finally runs into a singularity. At this stage it is necessary to go 
beyond the lowest order effective action for understanding how string theory cures the 
big-bang singularity, and for understanding the matching of pre-big-bang cosmology 
with standard Friedmann-Robertson- Walker (FRW) post-big-bang cosmology (the so- 
called 'graceful exit' problem [4-12]). One can imagine two possible scenarios. The first 
is that pertubative corrections succeed in turning the regime of pre-big-bang accelerated 
expansion into a decelerated expansion, and at the same time the dilaton is attracted by 
the minimum of its non-perturbative potential. This should take place before entering 
into a full strong coupling regime, so that perturbative results can still be trusted. 

In the second scenario the evolution proceeds toward the full strongly coupled regime. 
In this case one must take into account that at strong coupling and large curvature new 
light states appear and then the approach based on the effective supergravity action 
plus string corrections breaks down. The light modes are now different and one must 



1 



turn to a new effective action, written in terms of the new relevant degrees of freedom. 



In particular, at strong coupling D-branes [13] are expected to play an important role, 
since their mass scales like the inverse of the string coupling, ~ 1/g and they are 
copiously produced by gravitational fields ]TJJ (see also refs. [T7], |TH| for related 

approaches to singularities in string theory). 

In this paper we investigate how far one can go within the first scenario, analyzing a 
variety of perturbative string loop corrections. The motivation comes in part from the 
works |I| [TT|] , where the authors discuss general properties that the loop corrections 
should have in order to trigger a successful graceful exit. In particular, they find 
that the corrections should have the appropriate sign, so to induce violation of the null 
energy condition (NEC), but there should also be a mechanism that at some stage turns 
them off, otherwise the continued violation of the NEC produces an unbounded growth 
in the curvature and dilaton. In ref. [Tl[ they presented a model with loop corrections 
chosen ad hoc, both in the sign and in their functional form, just for the purpose of 



verifying explicitly on a toy model the general results of ref. JTU| , and they found that 
with appropriate choices a complete exit transition is indeed obtained. A transition 
to a regime of decreasing curvature was also obtained, with a potential chosen ad hoc, 
in the second paper of ref. [|l| . This shows that string loops can in principle trigger a 
complete graceful exit, but leaves open the question of whether this actually happens 
for the corrections derived from at least some specific compactification of string theory. 



In particular, the results of refs. [10, 11] show first of all that the corrections should 
have the appropriate sign, and this already is an interesting point to check against real 
string derived corrections, and furthermore they should also have a rather non-trivial 
functional form that suppresses them at strong coupling. 

String loop corrections have been much studied in the literature, especially for Z N 
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orbifold compactifications of the heterotic string [19-34], and we can therefore ask 
whether, at least in some compactification scheme, they fulfill the non-trivial properties 
needed for a graceful exit. In particular, corrections to the Kahler potential are known 
at all loops; this will be very important for our analysis, since in order to follow the 
cosmological evolution into the strong coupling regime, a knowledge of the first few 
terms of the perturbative expansion is not really sufficient, and one must have at least 
some glimpse into the structure at all loops. 

After discussing how far one can go within a perturbative approach, we will be able 
to shed some light on the question of whether perturbative string theory is an adequate 
tool for discussing string cosmology close to the big-bang singularity, or whether instead 
non-perturbative string physics plays a crucial role, as in the scenario discussed in 
ref. @. 

The paper is organized as follows. In sect. 2 we discuss the effective action of the 
heterotic string with a' and loop corrections for a orbifold compactification. In 
sect. 3 we use these corrections to study the cosmological evolution; we will introduce 
various corrections one at the time, to understand better their role, and we will find 
solutions that, in the string frame, smoothly interpolate between the pre-big-bang 
phase and a phase of decelerated expansion. In sect. 4 we examine the limit of validity 
of the perturbative approach. We point out that loop corrections induce an instability 
in the kinetic term of the dilaton, that vanishes at a critical coupling. Similarly to what 
happens in the Seiberg-Witten model, beyond this value of the coupling the effective 
action is more appropriately written in terms of variables related to the original ones by 
S-duality. In string theory this means that there is the onset of a new regime where the 
light modes of the effective action are given by D-branes rather than by the massless 
modes of fundamental strings. Sect. 5 contains the conclusions. 
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2 The effective action with loop corrections 



We consider the effective action of heterotic string theory compactified to four dimen- 
sions on a Zn orbifold, so that one supersymmetry is left in four dimensions, and we 
restrict to the graviton-dilaton-moduli sector. In a generic orbifold compactification 
there are the (1, 2) untwisted moduli fields C/j and the diagonal (1, 1) untwisted moduli 
fields Tj (non-diagonal moduli are included in the matter fields). The moduli fields Ui, 
determine the complex structure, i.e. the 'shape' of the compact space (they generalize 
the standard modular parameter of the torus); in various Zn orbifolds (see e.g. table 1 



in ref. p2|) the Hodge number = and therefore this shape is fixed and there are 
no moduli fields U. There are instead 3 diagonal moduli fields T,. We will neglect the 
fields U and we restrict to a common diagonal modulus, T« = T; it determines the 
overall volume of compact space. 

At the fundamental level string theory compactified on orbifolds is invariant under 
T-duality, which includes SL(2, Z) transformations of the common modulus T 

icT + d' {) 

with a,b,c,d G Z and ad — be = 1. These modular transformations are good quantum 
symmetries, and therefore they must be exact symmetries also at the level of the loop- 
corrected low-energy effective action. While at tree level the dilaton is inert under 
modular transformations, at one-loop the cancellation of the modular anomaly requires 
that the dilaton transforms as |19| 



S -> S + 2/dog(icT + d) , (2) 

where Re S = and is the dilaton field; k is a positive constant of order one which 
depends on the coefficient of the anomaly, see below. Then one easily verifies that 
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S + S + 2n\og(T + T) is modular invariant. We can therefore introduce a one-loop 
corrected modular invariant coupling g\ from [IS, 2C 

1 1 



, - n (S + S) + K\og(T + T) = e~* + Ka, (3) 



9l 2 

where we have defined the field a from ReT = (1/2)6°" (the factor 1/2 is not conven- 
tional but we found it convenient). Loop corrections to the effective action can be 



computed directly as an expansion in terms of the modular invariant coupling g\ []20 
As we see from eq. (||), 

and therefore an expansion in g\ provides a resummation and a reorganization of the 
expansion in e®. Note in particular that even when is large, the expansion in g$ is 
still under control if kct ^> 1, i.e. if nlogV ^> 1, where V is the volume of compact 
space in string units. 

We include in the action terms with two derivatives and terms with four derivatives, 
i.e. 0(a') corrections to the leading term. For both the two- and four-derivatives 
terms we include modular invariant loop corrections. We discuss separately the two- 
and four-derivatives terms in subsections 2.1 and 2.2. 

2.1 Terms with two derivatives 

In the Einstein frame, where the gravitational term has the canonical Einstein-Hilbert 
form, the action for the metric-dilaton-modulus system compactified to four dimensions 
is 

1 r M 1 

(5) 



S* = -( <*V=0 \\R ~ Kjd^Pzj 
a J iz 



Here z % = (S,T), K{ = d 2 K / dz l dzj and K is the Kahler potential. We use the 
conventions 77^ = (-, +, +, +), W vpa = d p Y^ a + . . . and R ua = R p upa . The superscript 
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E reminds that this action is written in the Einstein frame. The Einstein-Hilbert 
term ^—gR is not renormalized at one-loop |^6| and this non-renormalization theorem 
persists to all orders in perturbation theory around any heterotic string ground state 



with at least N — 1 space-time supersymmetry [34]. At tree level the Kahler potential 
is K tvcc = — log(5 , +5) — 3 log(T+T). It does renormalize, and at one loop, for heterotic 
string compactified on a orbifold, becomes [ITS 



# i-ioop = - log (S + S + 2k log(T + f )) - 3 log(r + f ) , (6) 

where k = 35 gs /(8tt 2 ). For instance for Z 3 orbifolds 5 GS = C(E 8 )/2 = 15, where C(E 8 ) 
is the quadratic Casimir of E 8 , and therefore k ~ 0.57. Eq. (Qj holds at one-loop, i.e. at 
first order in an expansion in 1/(S + S). However, the all-orders resummation implicit 
in eq. (@) is dictated by the fact that under the modular transformations given by 
eqs. (|1]J2|) the combination S + S + 2nlog(T + T) is invariant, and therefore respects the 
T-duality symmetry of the underlying string theory. In terms of the modular invariant 
coupling #q defined in eq. (§) we can write eq. (|) as 

^i-ioop = log (^pj - 3 log(r + f ) . (7) 

Eq. (0) is the leading term of an expansion in of the all-order Kahler potential. 
Indeed, the Kahler potential has been computed at all perturbative orders in ref. ||20|| , 
under the assumption that no dilaton dependent corrections other than the anomaly 
term are generated in perturbation theory. One defines implicitly the all-loop corrected 
coupling g 2 from 



1 1 2k , . 1 
— = — + -=-log(- 
9 2 9o 3 g 



2 - 2 + ^-log(-) + const. . (8) 



The all-loop corrected Kahler potential then reads |f20 

K = \o 



92 (l + ^ ' 



2 



31og(T + T), (9) 
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and at g\ <C 1 it reduces to eq. (0) plus terms 0{g\ log g^). 

The coupling g 2 is just the effective gauge coupling that multiplies the term F 2 after 



taking into account loop corrections |20| and, as we will see in the next section, it also 
multiplies the four derivative term. So the dilaton enters the action only through g 2 . 
We therefore define a new field <p from 

g 2 = (10) 

and we will treat it as our fundamental dilaton field. Therefore our loop-corrected 
two- derivative action is given by eq. (|5|) where now Zi = (S',T), Re S' = = g~ 2 
and K is given by eq. (H). 

In the following, it will be convenient to work in the string frame. In four dimensions, 
we define the string frame metric g^ v in terms of the metric in the Einstein frame g^ v 
from g^ v = g^e"^ (in four dimensions the dilaton is the same in the two frames). Note 
that we use ip rather than to transform between the two frames. At lowest order in 

of course this reduces to the standard definition, but beyond one-loop the definition 
in terms of ip is more convenient. 

Writing explicitly the kinetic terms of the dilaton and modulus field, (and omitting 
the superscript S for string frame quantities) the loop-corrected two- derivative action 
in the string frame reads 



2a' 

where we have defined 



R+(l + e^G(cp)) dptp&ip - ^aVa 



(11) 



2.2 Terms with four derivatives 
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2.2.1 The four-derivatives term at tree level 



For the four derivative term we find convenient to work directly in the string frame. 
At tree-level it can be written as [|35|, 



(<S0t.ee = — J d 4 x^=g - [R^R^ + bR^W + cR 2 ] (13) 

and, for the heterotic string, k = 1/2 (still, we display k explicitly in our equations to 
make easier the comparison with ref. ||; note however that we use the opposite metric 
signature compared to ref. M). The effective action can be obtained requiring that it 
reproduces the string amplitudes. This procedure fixes the coefficient of R ilvpa R^ vpa \ 
however, the coefficients b, c cannot be determined from the comparison with (on-shell) 
string amplitudes; this can be understood, for instance, showing that the coefficients 
b, c do not enter in the computation of amplitudes with three on-shell gravitons, while 
in the computations of four-graviton amplitudes they cancel between the contact and 



exchange graphs |31| 



A related source of ambiguity appears when one truncates the perturbative expansion 
in a' at any finite order. In fact, suppose that we choose somehow a value for b, c (two 
natural possibilities are either b = c = 0orb = —4:,c = 1, which forms the Gauss- 
Bonnet combination). Still, we can always perform a field redefinition that mixes 
different orders in a', e.g. g pv — > g fMU + a'(a 1 R pv + a 2 d^(f)d v (f)+ . . .), and (f> — > <f) + a'biR + 
. . .. This would not change the physics if one would be able to include all orders in the 
a' expansion, but it does make a difference if we truncate at a finite order in a'. For 
instance, — retaining only two- and four-derivative terms — this redefinition, applied to 
the two-derivative terms, generates new four derivative term, while the four- derivatives 
terms generates six- derivative terms which are truncated if we work at order a'. The 
terms which are generated at the four derivative level by a generic field redefinitions are 
just the terms R^R^, R 2 , plus other operators of the same order involving the dilaton, 
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like (<90) 4 , R^d^dyip, etc. This is another way to understand why the coeffcients of 
these terms cannot be fixed by the computation of string amplitudes. The term R^ upa 
is instead unaffected by field redefinitions, consistently with the fact that its coefficient 
is fixed by the comparison with string amplitudes. The study of the cosmological 
evolution using any specific action truncated at order a' should then be considered as 
only indicative of the possible cosmological behaviours. With suitable choices, one finds 
that the lowest-order pre-big-bang solutions are indeed regularized by the inclusion of 
a' corrections, and instead of running into the singularity, they are now matched (in 
the string frame) to a phase of De Sitter expansion with linearly growing dilaton |§. 
The effect on the cosmological evolution of these ambiguities have been studied in 
refs. ||, H 0. Here will take the point of view that, independently of these ambiguities, 
a solution that, thanks to suitably chosen a' corrections, approaches asymptotically a 
De Sitter phase with linear dilaton is a simple way to model a regularizing effect, which 
can have a different and deeper physical motivation: for instance, in ref. [0 we found 
that a similar transition to a De Sitter phase with linearly growing dilaton is driven 
by the formation of a gravitino-dilatino condensate. This mechanism is independent of 
the ambiguities discussed, and only depends on the dynamical assumption that such a 
condensate forms. 

Our choice for the form of the tree-level-four derivative term is the same used in 



refs. H. 11 



Octree = ^- ( ^ ] /,/', V ^,-' nf ils - (On) 1 • (11) 



2a' V 4 / 

In ref. || this form was obtained setting b = c = in eq. QT3D and then performing a 
field redefinition that generates R 2 and R 2 with the right coefficients to produce the 
Gauss-Bonnet combination R 2 GB = R\ vp <j~ + R 2 , and this also generates the term 
(<90) 4 . 
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Actually, our two-derivative action differs from that used in || already at tree level, 
because we have included the field a. The field redefinition that generates the Gauss- 
Bonnet term will also generate a number of four-derivative terms which depends also 
on da. We will neglect these terms because, on the one hand they make the action 
much more complicated, and on the other hand they are basically irrelevant to the 
dynamics, as will be clear from the results of sect. 3 and as we have checked on some 
examples. Actually, because of the ambiguities intrinsic in a truncation at finite order 
in a', it is not very meaningful to insist on any specific form of the action, and it is more 
important to look for properties shared at least by a large class of actions compatible 
with string theory. 

2.2.2 The loop-corrected four-derivatives term 

Let us first recall what happens in the slightly simpler case of a gauge coupling g a , 
where the index a refers to the gauge group under consideration. The coupling \j g 2 a is 
identified as the coefficient of (1/ '4)F" . At tree level, this term only appears when 
we expand in components the superfield expression (— 1/2) f a (S)W a W a , where W a is 
the chiral superfield containing F" ; fa{S) tree = S is independent of the gauge group, 
and \jg\ = ReS. At one-loop f a {S) gets a moduli-dependent renormalization [19-34] 

f a (S)^ loop = S + A a (T,T), (15) 

For orbifolds with no iV = 2 subsector, such as Z 3 and Z 7 , A a (T, T) = 5 a is a moduli- 
independent constant, and there is no moduli-dependent one-loop correction. Beyond 
one-loop, f a (S) is protected by a non-renormalization theorem |39| . 



Furthermore, at one loop a contribution to F^F^ comes from the anomaly: in 
fact, since the fermions in the supergravity-matter action are chiral, the tree level 
effective action leads, through triangle graphs, to one-loop anomalies. The type of 
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anomaly depends on the connections attached to the vertices. In particular, because 
of Kahler symmetry, chiral fermions are coupled to a Z/r-(1) Kahler connection, which 
is a non-propagating composite field. Modular transformations are a subset of Kahler 
transformations and therefore, if the theory is Kahler invariant, it is also invariant 
under modular transformations. Considering a triangle graph with one Kahler con- 
nection and two gauge bosons attached at the vertices, we get a mixed Kahler-gauge 
anomaly, proportional to F® u F aflu . The anomaly can be represented in the effective 
theory with a non-local term whose (local) variation reproduces the anomaly. Because 
of supersymmetry, this effective non-local term, when expanded in component fields, 
together with F®F a,J ' u also contains a term proportional to F^pa^v w j 1 j c j lj restricting 



to a common modulus, reads 119 



d 2 K 

K(T,T) - 2D" 1 — =W3„T 
v > J dTdT n 



\f; v F^. (16) 



with K(T,T) = — 31og(T + T)|] In the limit of constant T this term becomes local, 
and gives an additional moduli-dependent one-loop contribution to the gauge coupling 
g 2 a . Therefore, specializing for the moment to a Z 3 or Z 7 orbifold, 

]_ = s+i + T+ ^ 

Ta 2 

The variation of the term n\og(T + T) under modular transformation is just the 
anomaly, and the requirement of anomaly cancellation imposes the transformation 
law, eq. (|2|), on the dilaton field. Note that in the case of Z 3 and Z 7 orbifolds the 
moduli-dependent part of A a vanishes, and therefore it cannot contribute to the can- 
cellation of the anomaly; the cancellation comes entirely from the variation of S; so, in 

this case the anomaly must be independent of the gauge group, as indeed checked in 

1 A direct derivation of the term ~ F^F tJlV from Feynman graphs is quite subtle since, in the 
formulation of supergravity with the canonically normalized Einstein term, bosonic currents do not 
couple to the Kahler connection and therefore bosons cannot run into the triangle graph; this apparent 
puzzle is discussed in ref. M . 
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ref. [19]. Using the linear multiplet formalism [|D|, one realizes that this cancellation 
mechanism is just a four- dimensional version of the Green-Schwarz mechanism. 

At all loops, g 2 in eq. ( |17D is replaced by definition with the all-loop corrected effective 
coupling, which is the quantity that appears in the all-loop corrected Kahler potential, 
eq. © 0. 

Let us now discuss four-derivative gravitational couplings, which are the ones relevant 
for our analisys. In this case we have three couplings, multiplying R 2 , R 2 V ,R 2 . 
The one-loop renormalization of these coupling has been studied in refs. [26-31]. The 
situation is similar to the case of gauge couplings, and the contributions again come 
from the threshold corrections and from the anomaly (in this mixed Kahler 

-gravitational anomaly, i.e. a triangle graph with one Kahler connection and two spin 
connections attached to the vertices). There are however some complications: first 
of all, only one combination of these operators, corresponding to the the square of 
the Weyl tensor, i.e. to R\ vpc — 2R^ u R fll/ + (l/3)i? 2 , is obtained from a holomorfic 
function [[3(j] , and the other two independent combinations are not protected by a non- 
renormalization theorem. Since, in terms of superfields, the other two combinations 
that can be formed depend only on R 2 ^ u and R 2 , this means that R 2 lVfKr or R^ B are 
protected, while naked R^, R 2 terms are not. Furthermore, the ambiguity due to field 
redefinitions and truncation at order a' that we have discussed at tree level, persists 
of course at one and higher loops, so that terms like (dip) 4 , R^ u d^ipd u ip, etc., as well 
as term that depend on G(ip), can be generated with field redefinitions. It is then 
clear that the most general action is very complicated. We have chosen to focus on 
the loop correction to the same operator that we have considered at tree level, i.e. to 
the combination Rq B — (dip) 4 , neglecting naked R 2 ^, R 2 , (dip) 4 terms, as well as other 
terms that can be generated by field redefinitions. It is of course possible to extend 
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our analysis including other operators, but we believe that our choice is sufficient to 
illustrate the general role of loop corrections, while at the same time the action retains 
a sufficiently simple form, and in particular the equations of motion remain of second 
order. Our four-derivative action is therefore Si + «S n i, where 

1 (ka' s 



Sx = ^7 (x / ( 6 ~" + A((T) ) " {9ip) 



and iS n i is the non-local contribution from the anomaly, 

«- = 1 pr) tI" 1 ^^- 1 (Hf 9 "^) ■ < 19 > 

In the following we will neglect the non-local term. However, an effect of the anomaly 
is still present, because it has also produced the local contribution necessary to turn 
into the modular-invariant combination e _ip . Threshold corrections produce the 
function A, 

A(r, f) = log [(T + f ) |t7(«T)| 4 ] + 6 gI . (20) 
The constant 5 gT depends on the orbifold considered, and typical values are estimated 



in ref. [32 1. The constant b gT is related to the number of chiral, vector, and spin-| 



massless super-multiplets, Ns, Ny, N 3 / 2 respectively, by |$8| 

6 gr = ^(-3AV + N s ) - y (-3 + N 3/2 ) (21) 

and vanishes for orbifolds with no N = 2 subsector as Z3 and Z7] rj(iT) is the Dedeknid 
eta function, 

00 

V(r) = q 1/12 II (l - <l 2n ) > 1 = ^ ■ ( 22 ) 

n=l 

The anomaly produces also a term ~ R^ upiJ R pupa . However, below we will specialize 
to a metric of the FRW form, and in this background R ilvpa R! lVfXT vanishes identically 
(which allows us to look for solutions of the equations of motion with Im S = 0, 
Im T = H). 
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3 The cosmological evolution 

We now restrict to an isotropic FRW metric with scale factor a(t) = and Hubble 
parameter H = a/a = $. We use H to denote the Hubble parameter in the string 
frame. Another useful quantity is the Hubble parameter in the Einstein frame, He, 
related to H and ip by 

He = e^ 2 (H - \v) ■ (23) 

(Recall that we use (p rather than to move from the Einstein to the string frame). 
The shifted dilaton <p is defined by (p = ip — 3H. In the numerical analysis we will use 
units ka' = 1. 

At lowest order in both the a' and the loop expansion the solutions of the effective 
action come in pairs. For a = 0, they satisfy ip — (3 ± y3)H, or <p = ±\/3H, referred 
to as (±) branches respectively |]. The (+) branch, when H > 0, describes a Universe 
that starts from the low curvature regime and follows an accelerated superinflationary 
expansion. This is called a pre-big-bang type solution, and is characterized by <p > 
0. The (— ) branch, when H > 0, describes instead a post-big-bang evolution with 
decelerated expansion and (p < 0. In the present era we have a FRW decelerated 
expansion with stabilized dilaton, ip = 0, H > 0, and therefore (p < 0. A graceful 
exit is then obtained if the cosmological evolution connects the (+) branch with a 
(— ) branch solution. Necessary conditions for this to happen have been discussed 
by Brustein and Madden |TIJ, and it turns out that a rather non trivial interplay of 



different types of corrections is required. First of all, of course, we must move from 
the region of parameter space where (p > to the region with (p < 0. This can be 
accomplished with a' corrections |§. Furthermore, it can be shown that on the (+) 
branch He < (in fact, without corrections, He is monotonically decreasing and runs 
toward — oo at the singularity) and it is still negative at the branch change, i.e., when 
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(p = 0. To conclude the exit successfully, it is shown in ref. [10| that the evolution 
must proceed to a region where He is positive. This is in general non-trivial since at 
the beginning of the evolution He is negative and decreasing, and therefore a bounce 
in He is necessary. Thus, some new physical effect must be turn on to produce this 
bounce, and it is at this point that loop corrections are supposed to be crucial. It was 



found in ref. [l l| that indeed a one- loop correction with the appropriate sign turns 
the decrease of He into a growth and drives the solution into the region He > 0. 
However, this growth is unbounded, and the solution is now driven toward He — * +00. 
Therefore it is necessary that, after we have reached the region He > 0, some new 
mechanism turns on and kills the effect of the one-loop corrections. This was modelled 
in ref. [|ll|] multiplying the loop correction by a smoothed theta function. Then, after 
the evolution reaches the region He > and loops corrections are switched off, it 
becomes in principle possible to stabilize the dilaton with a potential. 

In this section we study the equations of motion, taking initial conditions of the 
pre-big-bang type; we will add various sources of corrections one at the time, in order 
to have some understanding of the role of the various terms, and we will compare with 
the above picture. 

3.1 The evolution without loop corrections 

First of all, we examine the behaviour of the system including a' corrections, but 
without the inclusion of loop corrections. In this case = (p and our action reads 

3 „ „„ ka' 



S = ^7 / d\^Tg e~* R + d^ip - \d,od»a + ^ (R 2 GB - (dip)' 



(24) 



If we neglect the modulus field a, this action reduces to that considered in ref. 
H and then we know that, starting from initial conditions of the pre-big-bang type, 
the solution has at first the usual pre-big-bang superinflationary evolution and then, 
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0.3 




(a) (b) 

Figure 1: (a) Evolution of H, ub for the tree- level system. Initial conditions at t — are: 
H(0) = 0.015, <p(0) = —30, <t(0) = 0.1; ip(0) is then fixed by the constraint equation, 
up = 0.16781 . . .. (b) The evolution of a. 



H 




9 

Figure 2: The evolution in the (H, (p) plane. The four lines, in order of increasing 
slope, are the (+) branch, the branch change line, the bounce line, and the (— ) branch, 
see the text. 



16 



when the curvature becomes of order one (in units ha! = 1), it feels the effect of the ol 
corrections. At this stage, instead of running into the singularity, it is attracted towards 
a fixed point with H and ip constants. This picture is not modified by the inclusion of 
a. In fact, writing also the equation of motion for a, one immediately sees that there 
is an algebraic solution of the equations of motion with b = 0, and H, (p constant and 
the same as in ref. [BJ, i.e. H = 0.616 . . . , <fi = 1.40 . . .. The numerical integration, see 
Figs, la, lb, shows that this solution is still an attractor of the pre-big-bang solution. 

For the discussion of the graceful exit, it is very convenient to display the solutions 
also in the (H, (p) plane, following ref. In this graph, shown in Fig. 2, four lines are 



of special interest. In order of increasing slope, the first line is the (+) branch of the 
lowest order solution (more precisely, this line corresponds to the lowest order solution 
only in the limit b = 0, and the deviation of the initial evolution from it that we see 
in Fig. 2 is due to a non- vanishing initial value of b). The second line corresponds to 
branch change, i.e. <p — 3H = 0. The third is the line where He = 0, and as found 
|T0f , it is necessary that the evolution crosses also this line to complete the exit. 



m 



Finally, we have the line representing a (— ) branch solution. We see from Fig. 2 that 
the lowest order solution ends up at a fixed point, after crossing the branch change 
line, but it is still in the region He < 0. 

The solution shown in this subsection can be considered as the starting point of 
our analysis; in the following subsections we will see how the various loop corrections 
modify this basic picture. 

3.2 The effect of the loop-corrected Kahler potential 

To begin our analysis we restrict to a Z3 orbifold, so that threshold corrections vanish, 
and we also neglect the non-local term. The action that we use in this section is 
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therefore 

S = — J d A x^~ge^ R + (1 + tfG(<p)) M ¥#V - -0„<70"o- + — (R 2 gb - (fy>) 



(25) 



with 



We again restrict to isotropic FRW metric and homogeneous fields and write the equa- 
tions of motion for the fields <p(t),a(t), /3(t). Taking the variation with respect to a, 
we get the equation of motion 

j t (e^a)=0. (27) 

Therefore, if we take as initial condition & = 0, o will stay constant. In this case the 
non-local term, eq. (|19| ), vanishes at all times, and therefore, for this initial condition, 
no approximation is made omitting it. 

Before starting with the full numerical integration it is useful to make contact with 



the general analisys of ref. JTT|. We therefore restrict to constant a (since a was not 
included in ref. [fLT| ) and we write the equation of motion obtained with a variation of 
the lapse function in the form 

QH 2 + ip 2 - QH(p = (p a , + Pq ) (28) 

where 

p a , = {ka')e-* (6#V - \^) (29) 

is the contribution of the a' corrections. The contribution of loop corrections is in the 
function p q which, from our action, turns out to be 

ft = -¥*?<?) = -¥> 2 (t)(3T^¥' (30) 
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In ref. [10| it was found that a graceful exit could be obtained with a loop correction 
that gives p q = —3f(ip)<p 4 , with /(</?) a smoothed theta function going to zero, for large 



<p, as e x%v . This form of the correction was just postulated in ref. [IT], but comparing 



it with the string result, eq. (BO), we find that, first of all, the sign comes out right, 
which is of course non-trivial. The dependence is ~ p> 2 rather than 4 since it comes 
from a correction to the kinetic term and, most importantly, its behaviour at large ip 
is different. In fact G(p) resembles a smoothed theta function, which is also a non- 
trivial and encouraging result, but it goes to zero only as e _</3 , which just compensate 
the factor e v in eq. (|28|) . We will see from the numerical analysis that this produces 
important differences compared to ref. [pj| . 

We now turn to the full numerical analysis, we restore a as a dynamical field and 
we set ka' = 1. The equations of motion obtained with a variation with respect to (p 
and (3 are, respectively, 

-6H(1 + H 2 ) + Cp(2 + 2e ip G + 3p 2 ) - 12H 2 - hs 2 - ^ - 6H* + 3Hp 3 + 
+6{l + e^G)Hp-p 2 {l -e v G') = 0, (31) 
4F(1 - Hip) - 20(1 + H 2 ) + 6H 2 - AHp + (1 - e^G)^ 2 - ^<p 4 - AH 3 <p + 
+2p 2 H 2 + hi 2 = , (32) 

and together with eq. ( BTT ) they determine the evolution of the system. The variation 
with respect to the lapse function produces a constraint of the initial data, 

6H 2 + <p 2 -mp-hy 2 = e^{p a/ + Pq ) (33) 

with p a i and p q given in eqs. (^,[3(]). The constraint is conserved by the dynamical 
equations of motion. We used this conservation as a check of the accuracy of the 
integration routine. Typically, the constraint is zero with an accuracy of 1CT 5 . 

The result of the numerical integration is shown in Figs. 3a, 3b. We see that at 
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(a) (b) 

Figure 3: (a) The evolution of 0, H including the all-order loop corrections to the 
Kahler potential; (b) the evolution of a. Initial conditions are the same as in the 
tree-level case and k = 0.57. 

first loop corrections are small and <p, H are the same as in Fig. la. At some stage 
loop corrections become important and the solution settles to a new fixed point, again 
with ip and H constant. Instead, at least on the scale used, the evolution of & is 
indistinguishable from the case without loop corrections, compare Figs. 3b and lb, 
because a is practically zero when loop corrections become effective. The change of 
regime takes place when the coupling g 2 is of order one, as can be seen from Fig. 4, 
where we expand the region in time where loop corrections become important and we 
plot H, ip and the coupling g 2 . 

From these plots, it might seem that after all the situation is not so different from the 
tree level evolution, because in both cases the solution in the string frame eventually 
approaches a De Sitter phase with linearly growing dilaton. An important difference 
however is found plotting the solution in the (H, (f) plane, see Fig. 5. We see in fact 
that the solution has crossed the line He = (and actually even the (— ) branch line) 
and therefore entered the region of parameter space where a graceful exit is in principle 
possible. Plotting the evolution of H E shows again that the loop corrections due to 
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Figure 4: if, and g 2 = as a function of time. Compared to Fig. 3a we have 
expanded the range of t where loop corrections become important. 

the Kahler potential produce a bounce in H E) see Fig. 6. 

Thus, loop corrections to the Kahler potential succeed in doing part of what loop 
corrections are expected to do, i.e. they produce a bounce in H E and move the solution 
into the region He > 0. However, we also want to obtain a solution with H, tp eventually 
decreasing and we want to connect this solution to the (— ) branch. We therefore turn 
to threshold correction to see if they can produce this effect. 

3.3 The effect of threshold corrections 

We now turn on the moduli-dependent threshold corrections, so that the action becomes 



S = — 




(34) 



The equations of motion are now (setting again ka' = 1) 



-6H(1 + H 2 ) + Cp[2 + 2e Lp G + 3(1 + A)0 2 ] - 12H 2 - ^& 2 - 
+ [3#(1 + A) + A]0 3 + 6(1 + e*G)H^ - 2 (1 - e^G') = , 




(35) 
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AH[1 - H(<p - A)] - 20(1 + H 2 ) + 6H' 
-AH\<p - A) + 2(0 2 + A)H 2 + h 2 = 



2 



mcp + (i - t*G)(p 2 



l + A 
4 








(36) 



and the constraint on the initial data is 




e^iPa' + Pq + Pga') , 



(37) 



where p a i , p q are given in eqs. (^,^) and 




(38) 



As initial conditions for a we take a value close to the self-dual point, <r(0) ~ cr sd = log 2, 
that is Re T ~ 1, and we take & small (consistently with the fact that the pre-big- 
bang evolution starts from the flat perturbative vacuum). We will discuss later the 
dependence on the initial conditions. With these choices, for a generic orbifold A(cr) 
turns out to be practically constant during the course of the evolution (and for a Z% 
or orbifold A(er) = 5 gT is exactly constant) and its value is determined by and 
S gr ; taking for instance 5 gT = 0, we have found nonsingular solutions in the range 
6 gr G [—20,0), which corresponds to A(a s d) G [—0.18,0). 

The evolution of the system under these conditions is shown in Figs. 7a, 7b. The 
behaviour of H, <fi is quite remarkable: threshold corrections turn the De Sitter phase 
with linearly growing dilaton into a phase with H, ip decreasing! At the same time 
the modulus a, and therefore the volume of internal space, shows a rather elaborate 
dynamics, see Fig. 7b. These figures refer to a Z e orbifold, for which k ~ 0.19. The 
same qualitative behaviour is obtained for a Z% orbifold, in which case A (a) = 5 gr is 
exactly constant, and the same results are also obtained for different, generic, values 
of K. 

The evolution in the (H, (p) plane is shown in Fig. 8, and we see that the solution 
approaches the (— ) branch. From this figure we also see that the solution approaches 
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Figure 7: (a) The evolution of H and (p with loop corrections to the Kahler potential 
and threshold correction switched on. The initial conditions are H(0) = 0.015, y?(0) = 
-30, o-(0) = 0.69, <r(0) = 0.001, and <p(0) = 0.07067. . . is then fixed by the constraint; 
the values of the parameters are n = 0.19, 6 gr ~ —4, 5 gr = 0. (b) The evolution of a. 
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Figure 8: The evolution in the (H, (p) plane. The straight lines are as in fig. 2. 

at first the tree-level fixed point discussed in sect. 3.1, then corrections to the Kahler 
potential and the threshold corrections become important about at the same time, so 
that after leaving this fixed point the solution deviates immediately from the behaviour 
that it has in the absence of thresholds corrections, shown in fig. 5, and it does not get 
close to the fixed point marked by a cross in fig. 5. Instead, if we do not include the 
corrections to the Kahler potential and we only switch on the threshold corrections, 
we found that the solution never crosses the bounce line H E = 0, and therefore the 
corrections to the Kahler potential are really an essential ingredient of our solution. 

Fig. 9 shows instead the evolution of the coupling g 2 , and we see that the curvature 
and the derivative of the dilaton start decreasing when g 2 ~ 1, so that when the solution 
is close to the (— ) branch we are already at large g, and at this stage non-perturbative 
effects are expected to become important. We will discuss this point further in sect. 4. 

Although it is appropriate to recall at this point that these results are obtained with 
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Figure 9: H, ip and g 2 against cosmic time. 

some specific choices of action and of initial conditions, it is certainly interesting to 
have at least an example of such a behaviour, with choices well motivated by string 
theory. To get some understanding of the dependence on the initial conditions we have 
run the integration routine for many different values of <r(0) and &(0). The shaded area 
in Fig. 10 is the region of the plane (a(0), <r(0)) where the behaviour is qualitatively 
the same as that shown above, while for initial conditions outside the shaded region 
the evolution in general runs into a singularity. Considering that a is at the exponent 
in Re T, the limitation on <r(0) is not particularly strong, while the required values of 
<t(0) are of the same order as the initial value of H . These initial conditions do not 
imply therefore any fine tuning. 

To have a better understanding of these solutions, it is also useful to display the 
corresponding Einstein-frame quantities. (We still plot them against string frame time 
t, but the same qualitative behaviour is obtained against Einstein frame time the 
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Figure 10: The shaded area indicates the region of initial conditions for which the 
system has a nonsingular evolution, and the dot corresponds to the value actually 
chosen in the solution displayed in Figs. 7-9. We have displayed only the part of the 
plane with o > cr s d = log 2, since modular invariance ensures that the figure is invariant 
under the transformation cr — > a sd — a. 

two are related by dt = (He exp(<^/2)). In Figs. 11a, lib we plot (pE = dip /(He and H E - 
The latter is particularly interesting and shows that in the Einstein frame our solution 
approaches asymptotically a De Sitter inflation. This is of course very different from 
the result of ref. || or of sect. 3.1, where De Sitter inflation takes place in the string 
frame. 
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Figure 11: (a) (pE and (b) H E against string time. 
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4 Transition to a D-brane dominated regime 



We now discuss the limitations on the validity of our solutions. As it is clear from 
Fig. 9, at large values of time we are deep into the strong coupling regime, g 2 ^> 1. 
Can we still believe our solutions? In our action we have included the corrections to the 
Kahler potential at all perturbative orders, while other operators, like R and -RfL^ are 
protected by non renormalization theorems. Therefore, despite the ambiguities that 
we have discussed for the four-derivative terms, due in particular to naked R 2 U and R 2 
terms, one might be tempted to argue that the solution is at least representative of the 
behaviour at strong coupling. However, this point of view is untenable, and at some 
point the perturbative approach itself breaks down. 

To understand this point, it is useful to work in the Einstein frame. The two- 
derivative part of our action then reads 



2a 

with 



S? = — d^x^g 



(39) 



Z v = l- 2e^G(<p) = l-2e*(—) / 6 + ^ 9 . (40) 
v Krj V 2 J (3 + nev) 2 v ; 

At weak coupling e v G(if) <C 1 and the kinetic term of the dilaton has the 'correct' 
sign. However, as e v — > oo, e v G((p) — > 3/2 and Z v < 0; Z v vanishes at a critical value 
of g 2 = given by 

At e v > g 2 it appears that the dilaton becomes ghost-like. The situation is quite similar 



to what happens in Seiberg-Witten theory [41], and the physical interpretation is the 



same. We can rescale the dilaton so that it has a canonically normalized kinetic term 
(—l/2)d ,J '<pd tl <p, and in terms of the rescaled dilaton the four-derivative interactions, 
and in general all interactions involving the dilaton, become strong as we approach 
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g c , and formally diverge at the critical point. This signals that the effective action 
approach that we have used breaks down and we must move to a new description, 
where the light degrees of freedom are different. In the Seiberg-Witten model the new 
weakly coupled modes are related to the original ones by an S-duality transformation. 
In string theory the S-dual variables are given by D-branes: for instance, in the Einstein 
frame, a fundamental string has a tension g 1 ^ 2 / (27ia') while a D-string has a tension 
g' 1 ' 2 / (2ira') [|T3 1 . This suggest that, if the cosmological evolution enters the regime 



e v > g 2 , the effective action approach that we have used breaks down, and we enter 
a new regime, which cannot be described in terms of a classical evolution of massless 
modes of a closed string; in this regime we must resort to a description in terms of 
D-branes. 

More precisely, the condition = identifies the critical point only if (p, H can be 
neglected. In fact, the equation of motion for (p in the Einstein frame reads (we insert 
for future use also a potential V(ip)) 

M v (p E = -3AH E ip E - V , (42) 

where 



1 - 2e^G(^) - 3A(a)0|J + . . . , (43) 
and 

A = 1 - 2e v G(<p) - A(a)<f 2 E + ... . (44) 

The dots denote tree-level a' corrections (which are neglegible at the later stage of the 
evolution). We recall that we found regular solutions for A < 0. So we see that, if we 
include the effect of the term |A|^|., the critical line is given by the condition = 
or 

Z v + 3\A\(p 2 E = 0, (45) 
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rather than = 0. Of course when ip^ approaches one we should at least include 
all higher powers in the a' expansion. More importantly, in the regime where H or ip 
approach one other restrictions on the validity of the effective action appear, and have 



been discussed in ref. embedding the 10-dimensional theory into 11- dimensional 
M-theory compactified on S 1 . In Fig. 12 (adapted from ref. [H]) on the vertical axis 
we show H, in the string frame. This is an indicator of the curvature and therefore 
of the typical energy scale of the solution. One might as well use 0, but of course 
precise numerical values here are not very important. In this graph we prefer to use 
the string frame quantity H because in this case the a' corrections become important 
when H ~ 1, while in terms of He this condition becomes e'^H^ ~ 1. 

The solid line H ~ 1/g separates the region where an effective 10-dimensional de- 
scription is possible, from the truly 11- dimensional regime. The region just above the 
line labelled lld-Sugra is described by 11-dimensional supergravity, while above the 
line labelled Dp-branes we are in the full M-theory regime. For a discussion of the 



crossover between these regimes we refer the reader to ref. [I3|. Of course, again, the 
position of the line separating the full M-theory regime from the 11-D supergravity 
regime is only indicative, and we have arbitrarily chosen its position so that it meets 
the curve H = 1/g exactly at H — g — 1. 

On the 10-dimensional side we have also drawn as a solid line the curve given by 
eq. (^), which is another critical line where a change of regime occurs. When ipE is 
not small, the form of this curve is only indicative. The position of g c depends on k, 
eq. fl4T|), and the graph refers to a Z e orbifold, k ~ 0.19. After the curve enters in the 
11-dimensional region it is not anymore meaningful. The label 'S-duality' means that, 
crossing this line, we enter a regime where the light degrees of freedom are related to 
the original ones by S-duality. On the same figure we display the solution of Fig. 7a, 
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Figure 12: The 'phase diagram' of M-theory compactified on S 1 . See the text for 
explanations of the various lines. The cosmological solution found in sect. 3 is marked 
by the arrows. 



labelled by the arrows. The solution for H will eventually decrease, but this only 
happens at very large values of g 2 (see Fig. 9), and we see that the solution enters the 
11-dimensional domain before it starts decreasing. At this point, it looses its validity. 

Finally, we found that it is not possible to stabilize the dilaton in our solution at 
a minimum of a potential, as could be generated for instance by gaugino condensa- 
tion [I3|. In fact at the later stage of the evolution the tree level a' corrections are 
neglegible, as we see in Fig. 7a, and M v ~ 1 — 2e ip G(ip)+3\ A\cp%. If we would stabilize ip 
around the minimum of the potential, it should first oscillate around the minimum and 
at the inversion points ipE = 0, so that here the coefficient of CpE in eq. ([42"D becomes 
~ Z<p. As shown in Fig. 13, this quantity is negative after we cross the He > line. 
As we discussed, this is not a problem for the consistency of the solution as long as (p 
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Figure 13: The evolution of Z v = 1 — 2e ip G((p) close to the point where He becomes 
positive, against string frame time. 

is not small (in fact, Fig. 12 shows that the limitation on the validity of the solution 
is rather given by the crossing into the 11-dimensional region), but it is clear that no 
consistent solution with ipE = can be obtained trying to stabilize the dilaton with a 
potential. In fact, if we try to force ipE to a small value, the coefficient of (pE in eq. ( f£2"D 
becomes approximately equal to Z v , which at this stage is negative. Therefore the evo- 
lution runs away from the minimum of the potential. Numerically, we have found that, 
including a potential in the numerical integration of the equations of motion, when the 
solution approaches the minimum of the potential the numerical precision, monitored 
by the constraint equation, degrades immediately and the solution explodes. 

Therefore, in our scenario, the problem of the dilaton stabilization can only be solved 
after the solution enters in the non-perturbative regime. 
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5 Conclusions 



In this paper we have tried to penetrate into the strong coupling regime of the cosmo- 
logical evolution derived from string theory. This regime is crucial for an understanding 
of the big-bang singularity in string theory, but since loop corrections do not tame the 
growth of the coupling while remaining within the weak coupling domain, it is clear 
that a knowledge limited to, say, one-loop corrections is of little use, and we really 
need to have at least a glimpse into the structure of the corrections at all perturbative 
orders. Luckily, for the effective action of orbifold compactifications of heterotic string 
theory, supersymmetry and modular invariance impose strong constraints on the form 
of the corrections at all orders. In particular, the kinetic terms of the dilaton is known 
exactly, while other operators, like R and R\ vpcn are protected by non-renormalization 
theorems. Therefore, in spite of some ambiguities in the choice of the four- derivative 
terms, present both at tree level and for their loop corrections, one can try to investi- 
gate string cosmology beyond the weak coupling domain, and to obtain at least some 
indications of what a well motivated stringy scenario looks like. 

As a first step, we have therefore tried to push this perturbative approach as far as 
possible, following the evolution even in the strong coupling domain g ^> 1. We have 
found solutions with interesting properties, that in the string frame start with a pre-big- 
bang superinflationary phase, go through a phase with H, (p approximately constant 
and of order one in string units, (a phase that replaces the big-bang singularity) and 
then match to a regime with H, p decreasing. Apart from their intrinsic interest, these 
solutions provide an illustration of the interplay between a' and loop corrections in 
string cosmology, and give an explicit realization of the general picture emerged from 
the works [1,5-12]. Probably the main element that is missing from this part of the 
analysis is the inclusion of non-local terms. These might model the backreaction due 
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to quantum particles production, which might play an important role in the graceful 
exit transition P2"| . Unfortunately, these are quite difficult to include in a numerical 
analysis. 

Despite some nice properties, the cosmological model that we have presented still 
have some unsatisfactory features, and in particular the dilaton could not be stabilized 
with a potential, and so this model cannot be the end of the story. 

On the other hand we have found that, if we look at our solution from the broader 
perspective of 11-dimensional theory, it ceases to be valid as soon as we enter into 
the strong coupling region, even if one includes perturbative corrections at all orders. 
Thus, we think that our analysis reveals quite clearly the direction that should be 
taken to make further progress. As already discussed in ref. [Q, when we move to- 
ward large curvatures we meet critical lines in the (H, g) plane, beyond which D-branes 
becomes the relevant degrees of freedom. Here we have found another critical line at 
strong coupling; beyond this line the light modes relevant for an effective action ap- 
proach are obtained by an S-duality transformation and are therefore again naturally 
interpreted as D-branes. The combination of these critical lines, shown in fig. 12, and 
the behaviour of our solutions, also displayed on the same graph, suggest that the 
evolution enters unavoidably the regime where new descriptions set in. The under- 
standing and the smoothing of the big-bang singularity therefore requires the use of 
truly non-perturbative string physics. 
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